US007817087B2 


(i2) United States Patent 

Cohen et al. 


(io) Patent No.: US 7,817,087 B2 

(45) Date of Patent: Oct. 19, 2010 


(54) METHOD AND APPARATUS FOR RELATIVE 
NAVIGATION USING REFLECTED GPS 
SIGNALS 

(75) Inventors: Ian R. Cohen, Eldersbuig, MD (US); 

Gregory J. Boegner, Jr., White Hall, 
MD (US) 


(58) Field of Classification Search 342/353, 

342/357.06, 357.08, 357.13, 357.15, 357.16; 
455/12.1, 13.1, 13.3; 701/213,215 
See application file for complete search history. 

(56) References Cited 

U.S. PATENT DOCUMENTS 


(73) Assignee: The United States of America as 

represented by the Administrator of 
the National Aeronautics and Space 
Administration, Washington, DC (US) 


6,653,972 B1 11/2003 

2005/0192719 Al * 9/2005 

2006/0077094 Al 4/2006 

2006/0082496 Al 4/2006 

* cited by examiner 


Krikorian et al. 

Sheikh et al 701/13 

Eisenstein 
Winternitz et al. 


( * ) Notice: Subject to any disclaimer, the term of this 

patent is extended or adjusted under 35 
U.S.C. 154(b) by 0 days. 

(21) Appl.No.: 12/116,518 

(22) Filed: May 7, 2008 

(65) Prior Publication Data 

US 2009/0278736 Al Nov. 12, 2009 

(51) Int.Cl. 

G01S 19/42 (2010.01) 

G01S 19/51 (2010.01) 

(52) U.S. Cl 342/357.25; 342/357.34 


Primary Examiner — Dao L Phan 

(74) Attorney, Agent, or Firm — Matthew F. Johnston 

(57) ABSTRACT 

A method and system to passively navigate an orbiting mov- 
ing body towards an orbiting target using reflected GPS sig- 
nals. A pair of antennas is employed to receive both direct 
signals from a plurality of GPS satellites and a second antenna 
to receive GPS signals reflected off an orbiting target. The 
direct and reflected signals are processed and compared to 
determine the relative distance and position of the orbiting 
moving body relative to the orbiting target. 
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METHOD AND APPARATUS FOR RELATIVE 
NAVIGATION USING REFLECTED GPS 
SIGNALS 

GOVERNMENT INTEREST 5 

The embodiments of the invention described herein were 
made by employees of the United States Government, and 
may be manufactured and used by or for the United States 
Government for governmental purposes without payment of 10 
any royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 15 

The present invention relates to navigating a moving body 

with the use of GPS signals and more particularly to the use of 
reflected GPS signals to passively navigate an orbiting body 
towards an orbiting target. 

2. Description of the Related Art 20 

Spacecraft require equipment that calculates the orbital 

position, velocity, time, and attitude. The advent of the Global 
Positioning System (“GPS”), increased the desirability of 
using onboard processing systems for spacecraft position 
determination. 25 

The GPS satellite system is a collection of satellites that 
can be used for missile, satellite, aircraft, and terrestrial navi- 
gation. Each GPS satellite broadcasts its own ephemeris and 
time, thereby allowing a GPS receiver to determine its posi- 
tion. Typically, a GPS receiver calculates its position from the 30 
simultaneous observation of any four GPS satellites in order 
to calculate its position. A civil grade GPS receiver can accu- 
rately determine its position within 20 meters, determine its 
velocity within 0.6 meters/second, and the current time with 
a 10-nanosecond accuracy. At least twenty-four valid GPS 35 
satellites are in the GPS constellation at any given time. 

A GPS receiver is an autonomous instrument that trans- 
forms signals from GPS satellites into point solutions for 
spacecraft navigation. Current GPS receivers have a radio 
frequency section for receiving and converting signals 40 
received from a spacecraft’s antennas. The digitized signals 
are then forwarded to one or more correlators, controlled by 
the receivers own processor. The correlators look for matches 
between the incoming signal and the C/A code for different 
satellites . When a satellite lock occurs, or the incoming signal 45 
matches an internally generated pseudo random noise (PRN) 
code, the receiver’s processor is notified. The processor con- 
tains executable code to generate a pseudo-range or a line-of- 
sight distance to the satellite. The processor may also contain 
the executable code of an orbit propagator. The pseudo-range 50 
is a measurement input to the navigation filter that calculates 
point solutions for determining the orbit of the spacecraft. 
While such GPS systems may be employed to position and 
track an orbiting moving body such as the Space Shuttle, such 
systems do not provide the ability to determine the position 55 
relative to another orbiting body. Conventional GPS receivers 
can not passively provide relative navigation or relative posi- 
tion data which are necessary during formation flying, dock- 
ing, or other approach of space craft to proximity to another 
orbiting body. 60 

The space environment poses additional obstacles to that of 
land based tracking systems. Because of the speeds associ- 
ated with spacecraft motion, GPS satellites are in view for 
much shorter time periods than when viewed from a ground 
reference. The appearance and disappearance of GPS satel- 65 
lites from view causes their output signals to slew through the 
entire range of the 45 -kilohertz Doppler shift. Typical space 


2 

craft formation flying and autonomous rendezvous missions 
rely on active transmissions schemes such as RADAR (RA- 
dio Detection and Ranging) and LIDAR (Light Detection and 
Ranging) which require specialized hardware requiring addi- 
tional mass and power consumption. With the extreme weight 
and power restraints associated with space flight, there is a 
particular need for a reliable navigation/positioning system 
which reduces weight and power consumption requirements, 
as well as taking advantage of free GPS signals. 

The use of Reflected GPS signals per se is also known. U.S. 
patent application publication 2006/0077094 entitled Infor- 
mation “Information Gathering Using reflected Signals” dis- 
closes a method of using reflected signals in a land based 
application to monitor the presence or absence of an object 
such as if a vehicle occupies a parking space and is hereby 
incorporated herein in its entirely. This system is not only ill 
suited for orbital applications, but does not make use of the 
reflected signal to navigate to a target. 

It is also known to use Bistatic Radar Systems where a 
radar transmitter located remotely from the projectile (such as 
onboard a ship) illuminates the target and the reflected returns 
are received by a receiver located on the proj ectile. The track- 
ing data from the radar measurements are then used to calcu- 
late the proper guidance signals to direct the projectile to the 
target. An example of a bistatic radar system is discloses ill 
U.S. Pat. No. 6,653,972 which in incorporated herein it is 
entirely. These systems require the use and control of an 
active radar transmitter and do not provide a suitable global 
positioning function. Rather, these systems are employed to 
simply acquire and hone in on a particular target. Therefore, 
such Bistatic radar systems are not suited for orbital position- 
ing and navigation or have the ability to rely on a passive 
autonomous system. 

One of the objects of the present invention is to utilize 
reflected GPS signals to provide passive autonomous relative 
navigation of an orbiting spacecraft towards an orbiting body. 

SUMMARY OF THE INVENTION 

The present invention is directed to a method and system to 
passively navigate an orbiting moving body towards an orbit- 
ing target using reflected GPS signals . A pair of receive anten- 
nae are employed. One to receive direct signals from a plu- 
rality of GPS satellites and a second antenna to receive GPS 
signals reflected off an orbiting target. The direct and 
reflected signals are processed and compared to determine the 
relative distance and position of the orbiting moving body 
relative to the orbiting target. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram depicting the relative geometry 
between an orbiting spacecraft, target and satellite. 

FIGS. 2-3 are diagrams depicting the earth centered inertial 
coordinate frame. 

FIG. 4 is a diagram depicting the visibility of GPS signals 
to a body in orbit. 

FIG. 5 is a block diagram of a system to receive and process 
reflected signals according to a preferred embodiment of the 
present invention. 
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FIG. 6 is a flow diagram depicting the algorithm to receive 
and process reflected signals according to a preferred 
embodiment of the present invention. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT 

A preferred embodiment of the invention and the various 
features and advantageous details thereof are explained more 
fully with reference to the accompanying drawings and 
detailed in the following description. It should be noted that 
the features illustrated in the drawings are not necessarily 
drawn to scale and descriptions of well-known components 
and processing techniques are omitted so as to not unneces- 
sarily obscure the embodiments of the invention. The 
examples used herein are intended merely to facilitate an 
understanding of ways in which the embodiments of the 
invention may be practiced and to further enable those of skill 
in the art to practice the embodiments of the invention. 
Accordingly, the examples should not be construed as limit- 
ing the scope of the embodiments of the invention. 

The present invention is directed to a method and apparatus 
for passively navigating an orbiting space craft towards and 
orbiting target using reflected GPS signal. FIG. 1 depicts the 
relative geometry between the space shuttle, Hubble space 
telescope and a GPS satellite. By comparing both signals 
received directed from the GPS satellite and signals reflected 
off the Hubble, a navigation system employed on the space 
shuttle can passively and autonomous determine its position 
relative to the Hubble thereby providing navigational aid 
passively and autonomously. Acquisition of and tracking of 
GPS signals are well known in the art and need not be dis- 
cussed in detail. U.S. published patent application US 2006- 
0082496 to Wintemitz et al. discloses a radiation-hardened 
fast acquisition/weak signal tracking system and method and 
is hereby incorporated herein by referenced in its entirety. 

The actual GPS measurements recorded from any of the 
constellation of satellites are the pseudorange and psuedor- 
ange rate. The pseudorange, p, is based on the time of flight 
between the GPS receiver and the GPS transmit antenna. This 
is used to form a range measurement, and thus solve for a 
position fix when enough satellites are available. The other 
measurement available is pseudorange rate that is essentially 
the rate of change along the line of sight vector to the satellite. 
The pseudorange rate, dotp, can also be viewed as the Dop- 
pler shift on the signal. The pseudorange rate equation is 
equivalent to the carrier rate, <E> which tends to be a less noisy 
measurement. 
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The GPS satellite clock error is the residual error from a 
second order polynomial fit, adjusting for relativistic effects, 
the clock bias, clock drift, and frequency drift. The Satellite 
Clock error is on the order of 3.0 m. The satellite ephemeris 
error is caused from the deviation in the current satellite 
ephemeride stored onboard the GPS satellite with the satel- 
lite’s true position. These ephemeris values are used to fonn 
the line of sight vectors from each SV to the user, thus each 
perturbation adds error into the pseudorange. This is normally 
on the order of 4.2 m for a typical ground based GPS receiver. 

The GPS system also experiences relativistic effects due to 
the gravitational potential differences. Space users also expe- 
rience ionospheric effects. This is generally seen as a bias in 
the measurement. The ionospheric effect is modeled off of the 
path length of the signal traveled through the ionosphere. 
There are also measurement errors from within the receiver. 
These are from the tracking loops internal to the GPS receiver. 
The dominant sources of error introduced here are the thermal 
noise jitter and dynamic stress error. The secondary sources 
are given by the hardware and software resolution, and clock 
drift. The last sources are multipath and shadowing effects. 
For the purposes of rendezvous, we are actually interested in 
tracking these multipath sources in order to solve for a relative 
position. The GPS pseudorange error budget is on the order of 
8.0 m for a standard GPS receiver located on the ground. 

The direct pseudorange rate is a measurement of the rate of 
change of the psuedorange between the user and the GPS 
satellite. It is given by 

P=(*>0’16s +e 

That is the rate of change of the psuedorange along the line 
of sight vector and the user. The error is mainly determined by 
the clock drift component. Generally these measurements are 
noisy compared to that of the pseudorange. 

The constructed measurements for the bistatic ranging 
problem are the differenced reflected psuedorange and the 
reflected psuedorange rate. The reflected psuedorange, is 
differenced with the direct signal psuedorange to remove the 
common mode errors from the measurement. Some of the 
common mode errors removed are front-end jitter, receiver 
bias, thermal noise, and line bias, leaving us primarily with 
ionospheric effects, multipath and jitter error contributions. 
The psuedorange rate is the additive Doppler shifts along the 
line of sight vectors between Hubble and the Shuttle, and 
between Hubble and the GPS constellation. 

This difference reflected pseudorange measurement is con- 
structed from the difference in the direct path signal and the 
reflected signal. The equation is given by: 


The measured pseudorange is the difference between the 
system time at which the signal left the GPS satellite and the 
system time at which the signal reached the user. This is given 
by 

p=r+c(t u +bt+bt D ) 

r is the geometric range is given by r=c(T M -TJ, which is the 
exact time of flight between the user and the GPS satellite. t u 
is the receiver offset from system time, b t is the offset of the 
satellite from system time, and h tD is due to other measure- 
ment errors. 

The GPS constellation is segmented into 3 parts, the user, 
the control, and space segment. The three main sources for the 
space segment are satellite clock stability, clock perturba- 
tions, and selective availability. For the control segment, it’s 
primarily ephemeris prediction error. For the user segment it 
is ionospheric delay, multipath, receiver noise, and resolu- 
tion. 


P pr P/? — Pi) 

Expressing as a geometric relation we arrive at: 

P P r=r sh +r hg -r gs +e 

That is the sum of the geometric range from the Shuttle to 
Hubble, Hubble to the GPS satellite, and the Shuttle to the 
GPS satellite, plus the residual errors. 

The psuedorange rate can be seen as a Doppler shift on the 
signal. The reflection can be viewed as a retransmission, 
making it an additive process. 

p r Hv s -v h )ids sh Hv h -v g )[ds hg 

Where the los^ is the line of sight vector between the 
Shuttle and Hubble, and los^ is the line of sight between 
Hubble and a GPS satellite. The primary sources of error 
contributing to the reflected psuedorange rate are the same of 
that as the direct measurement, however, the quality of the 
incoming signal is degraded due to the reflective process and 
the actual effect is unknown. 
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Typical GPS receivers use a simple least squares algorithm 
to determine position and time. It requires at least 4 visible 
satellites to solve for the four unknowns x,y,z, and the user 
clock offset from the GPS constellation. In addition to posi- 
tion estimates, some also do perform a finite difference of the 5 
positions to solve for velocity as well. 

In more advanced receivers, a Kalman filter may be imple- 
mented in order to utilize the dynamic information, as well as 
the measurements to improve the accuracy of the receiver. On 
orbit, the GPS receiver will experience large dynamic stress 
and motion with respect to the GPS satellites, it is almost 
necessary to filter this heavily to provide a good navigation 
solution. The method of the present invention implements an 
extended Kalman filter which will now be described. 

The Kalman filter represents an optimal filter for a linear 
system with additive gaussian noise with zero mean. Since the 1 5 
dynamics of an orbiting body are well understood and their 
equations are well formulated, it is possible to exploit that 
information when forming the Kalman filter when applied to 
the relative navigation problem. The canonical form of the 
Kalman filter is the optimal filter for a linear system, incor- 20 
porating a blending of known dynamics, with that of the 
measurements. The discrete Kalman filter can be broken up 
into two distinct steps, the time update and measurement 
update. 

The dynamical equations in a standard linear system are in 25 
the form: 

X=Ax+Bu+d 


y=Hx+n 

30 

These equations are continuous and need to be discretized 
in order to implement them in the Kalman filter. 

This discretization is given by: 

% + i=^+r*%+4 

35 

y k =Hx k +n k 

This relates timestep k-1, to timestep k. In order to calcu- 
late the state evolution between timesteps, for a linear system, 
the state transition matrix O is given by: 

40 

® k = e AAT 

Where 


AT-t k+l -t k . 

The evolution of T is a bit more complex, it is given by: 45 
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The error covariance propagation step is given as follows: 

P k+l -=<s>kP k ® k T +Q 

The first step of the measurement update is to form the 
Kalman Gain, K^. The Kalman gain is the optimal estimate to 
minimize the errors along the diagonal of the covariance 
matrix. 

K k =P k ~H k T (H^H k T + R)- 1 

The next step is to apply the gain to the measurement error. 
The H matrix is the measurement connection matrix, it con- 
nects the states with the measurement at time- step k. 

Xk= x k~+Kk(Zk-HkX k -) 

The final step is to correct the error covariance for time-step 
k. 


P k =(I-K k H)P k ~ 

While die equations are relatively straight forward in the 
discrete Kalman filter, they require a measurement model and 
the plant dynamics are linear systems. Unfortunately, in the 
real world most systems are nonlinear. In particular, GPS has 
a non-linear measurement model. 

The solution to a nonlinear measurement model can be 
done by approximating Y K . Given a nonlinear measurement 
model Y=h(x) 


y k = H k x k 


H k 


dh 

dx 


\x k 


This approximation is called an extended Kalman filter. 
The last two parts of a Kalman filter are the measurement 
covariance, and process noise matrices, R and Q respectively. 
The measurement covariance matrix consists of the variance 
of noise on the signal. From the canonical form: 

yk=HjPk+”k 


it can be shown that 
R=E[W] 

Assuming uncorrelated gaussian white noise with zero 
mean, R is given by: 


r T 


C pi nxn 

e AAT Bu(T)dr 

Jo 

R = 

Q*» (T 2 plran 


50 


This model only concerns the unforced dynamics, so with- 
out loss of generality, B is assigned 0 and T is dropped. Now 
these equations can be applied to the Kalman filter form. 

The first step in the Kalman filter is to propagate the mea- 55 
surement from time-step k to the time-step k+1 . This is done 
by multiplying by the state transition matrix, or integrating 
from 


The process noise matrix is given by the variance un- 
modeled forces in the dynamic model. Assuming we have a 
model: 

** =< £***+r *%+4 

noise actually enters the system as a disturbance on the input. 

Q=r^[dA r ]r/ 

Unfortunately, the equation to describe the matrix Q is 
difficult to obtain because it requires calculation of T. In the 
simulations, Q is treated as a free variable, adjusting it to give 
the optimal performance, and reduce the process covariance 
bounds such that they were roughly equal to the measured 
variance in R. 

In the case of relative navigation, it is an essential to be able 
to solve for your relative position as well as possible. One of 


60 

4 + r=^A+U% 

Propagating the error covariance matrix. 

P is defined as 

P=E[(x-x)(x-x) r ] 65 

it represents the variance in the estimator errors. 
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the best ways is to improve your solution is to have an accu- 
rate dynamic model. This way, one can rely on the dynamic 
model to propagate through measurement outages and noisy 
measurements. This is especially true in the case of applying 
reflected GPS measurements. 

Defining the necessary components to construct the filter is 
necessary before beginning formulation of the Kalman filters 
for the relative navigation problem. Given a dynamical sys- 
tem: 


8 

typically one spacecraft is chosen to be the reference orbit. 
This orbit is ideally circular. In this case, the orbit chosen is 
the space shuttle. This is denoted by r tg£ , the other spacecraft 
is typically called the interceptor or r zwr Hill’s equations, also 
5 known as the Clohessy Wiltshire equations which are the 
governing forces for relative two-body motion. These equa- 
tions are well known and have an explicit solution. 

x-2wy-3 m^jc+/x=0 
10 y+2wx+fy=0 


y=h{x) 

It is thus necessary to define a state vector, X. The state 
vectors used in an absolute model, are the absolute positions 
and velocities of the Shuttle and HST, as well as the bias state 
for the navigation solution. Typically a GPS receiver would 
include the bias drift as another component. For this simula- 
tion, since the rendezvous is over a short period of time, bias 
drift can be neglected. 

T 20 

x=[r s ,v s ,r h ,v h ,b] 

For a completely relative model, Hill’s equations are used 
to propagate the states. This requires formulating the problem 
entirely in a relative frame. The states in this case are: 

X^Avf 25 

In this formulation, A r is the relative position in hill’s frame 
and A v is the velocity. Typically the dynamic model is a 
nonlinear equation given by a f(dot)=f(x). In this analysis two 
different model are examined. The first model is based purely 30 
on the 2 body equation of motion given by Hill’s equations. 
This model is a linearized representation of the relative space- 
craft dynamics and has the form 

x=Ax. 

35 

where A is given by the partials of f(x) evaluated along a 
reference trajectory. 



Hill’s model is conveniently linear, simplifying the appli- 
cation of the Kalman filter. Hill’ s model is only valid for close 
proximity maneuvers and requires a special coordinate frame 
to operate in, centered on a circular orbit for one of the 
vehicles. 

Hill’s equations of motion are derived in the R, S, W 
coordinate frame as depicted in FIG. 2. 

R is the radial direction from the center of the earth, 

S is the projection of the velocity in the plane perpendicular 
to R, and W is the direction of angular momentum perpen- 
dicular to the orbital plane. 

The transformations from the IJK coordinate frame are as 
follows: 55 



S=WxR 

RxR 

w = 

\Rxk\ 


z+w 2 z+fz = 0 

or, by defining 

x r =(x,x), then provides the form x=f(x). Now, taking par- 
tials with respect to the states, arriving at the canonical linear 
system form x=Ax. 
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Hill’s equations being in a linear form fit well into the 
linear systems toolset. If we assume a constant A over each 
discrete time step we can approximate: 


= * (AAt) 

df I 

A = ~r\ 

dx |* 4 


This formulation provides a state transition matrix in linear 
form that can be easily applied to both the state propagation 
steps and the covariance propagation. Although the relative 
dynamics are minimal in comparison between spacecraft, 
their actual perturbations on the dynamics are quite signifi- 
cant in low earth orbit. In addition to increasing the fidelity of 
the dynamic model, it allows customization of the orbits and 
adds other effects such as thruster firing, and differences in 
their drag coefficients. 

The coordinate frame used in the high fidelity dynamic 
model is the Earth Centered Inertial frame. This reference 
frame was chosen because it greatly simplifies the equations 
of motion. The ECI coordinate system is also known as the 
IJK system. Referring to FIG. 3 the I axis is pointed towards 
the vernal equinox, the J axis is 90° to the east along the 
equatorial plane, and the K axis extends through the North 
Pole. 

The lull nonlinear dynamic equations are given by: 



1 A 

-C D -pV(x + u) E y) 

2 m 

1 A 

~C D -pV(y + a) E x) 

2 m 


1 A 
C D -pVz 

2 m 


Where r and r are the position and velocity vectors in the 
IJK reference frame. Since Hill’s frame is a relative frame, 


65 


The J 2 term is the zonal harmonic constant given by 
1 082.7x1 0 -6 . co £ is the rotation rate of the earth, R E is the 
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radius of the earth, m is the objects mass, A is the cross 
sectional area, and V is the relative wind. The relative wind is 
given by the equation 


10 


dr y 

dx 


= 3^|l-l 
2 


-continued 
1 + 




V = + 2oj Er x r y + r 2 y - 2 u) E r y r x + w|r 2 + r\ 


The value of p is given by an exponential model derived in: 


r=v and v=g(r,v) 

In this formulation, g is the nonlinear force model for the 


orbital dynamics. Now, let the state vector x= 

x=[v,g{r,v)] T 


|r,v| 


1 A w x(y + u E x) l /i 

- C D —pv - -C D —p 

2 m rH 2 m 


Voj e + 


u E {y + ^Ey) 
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Where p 0 =3.275e-12 kg/m 3 , a is the current altitude, a 0 is 
the reference altitude, 400000 m and H is the scaling height 
58515m. Using the same from the simple GPS Kalman filter, 
choosing r as position and v as velocity, the equations it is now 1 5 
appropriate to write these equations in terms of first order 
derivatives, to fit the canonical form, are 


dr y 

dy 


3 /Re\ 2 ( z 2 

1 ~At) ’t*- 1 


+ 3 


py 


5 / R E \ 2 / z} 


20 


1 A w y(y + a> E x) 1 A oj E {y + oj e x){x- cj E y 

2 Cd m^ V ~rH + 2 C °m P V 


dr y 1 A (y - io E x)(x - io E y) 

= -~C D —p— 
dx 2 m 


<- o E y ) 


dr y 1 A . 

jf=-2 c °A v+ 


V 

(y + (x) E xY 
V 


dr y 1 A Z(jf + o) E z) 

TZ = ~n C D~P 77 

dz 2 m V 




r e 


3 + 


f(x) = 


1 A Zx 1 A (o E z(y + o) E x ) 

~C D —pV — - -C D —p 

2 m rH 2 m V 


30 


35 


dr z 

dz 


dr z 

dz 


.3«*|l_!^| |7^-3 




1 „ A zy 1 A ( 0) E z(x - 0) E y ) ) 

~C D -pV— --C D -p\ 

2 m rH 2 m \ V ) 


4li-3*(^ff4-3|| + 


The equations are in the appropriate form to apply the 
Kalman Filter. Since the Kalman filter’s covariance propaga- 
tion is a linear process, it is necessary to formulate the partial 
derivatives with respect to the state vector. Although the par- 
tials are not used to propagate the dynamics, they are needed 40 
when forming the state transition matrix for the covariance 
updates. 


3^|l-U 2 (^) 2 f 7 4- 5 )] + ic D V^ 


dr z 1/4 Z(x + co E y) 

= ~n^D—P 7} 

ok 2 m V 

dr z 1/4 Z(y + oj e x) 

~7T = — o '-'D P 7} 

dy 2 m V 

dr z 1 A ( 

-zf = -o C D-P\ V+ 17 
at 2 m y V J 


dr x = 1 dry = i dr z 
dx dy dz 


45 


dr x 

dx 


dr x _~P*y 

dy r 5 


1 A x{x + u) E y) 1 A x{x + io E y){y - u) E x) 50 

-C D -pv — + - C D —pV 

2 m rH 2 m V 




2 '4 


, 1 A y{x + oj E y) 
+ 7>CD-pV- 


rH 


1 A 
-C D —p 

2 m 


Vu> E 


w E {x + (o E yr 


dr x _ „PXZ 
d Z ~ I s 




1 A 

77 C D-PV 


V 

Z(x + (x) E y ) 

rH 


dr x 1 A 

T7 = --xCd—P 

dx 2 m 




{x + oj E yf 


dr x 1 A {y - oj e x){x - oj E y) 

~n~T - — 77 '-'D P 7} 

dy 2 m V 


dr ■ 


1/4 z{k + oj E y) 

~7T = ~ n I'D P 77 

dz 2 m V 


Hill’s equations of motions are conveniently in linear form: 
X(^)=A x +B m . The state transition matrix can simply be 
approximated by the matrix exponential 


Where A is the linear system A matrix, and T is the time 
interval. In this particular implementation the time step is on 
the order of 1 second, thus T=1 This is both used in the 
dynamics, and the covariance propagation. In the full nonlin- 
55 ear case estimates would diverge after a few iterations using a 
first order approximation of the system. Thus the state is 
integrated over time steps of one second. The covariance 
matrices are propagated linearly, using, 

Pk+i =< &(h+i>tk)P k®(tk+i A) T+ Q- 
Solving for O in a given linear system, we have: 


60 


65 


®(t k+l ,t k H tk tk+lA A)®dx 

and arrive at an approximate 
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Since A is not constant over the time interval t k+1 , t*, we 
cannot rely on the matrix exponential approximation for this 
integration. 

In order to incorporate the variation of the states between 
time steps, the added dynamic equation <b=Ad> is included 5 
within the dynamic model propagation. This allows integrat- 
ing the dynamics and solving for (|) between time steps, while 
varying A with the current state estimates . In order to solve for 
(|), the initial condition vector of the states and the identity 
reshaped into a vector is required. The resultant (|) at the end of 10 
the integration is then used in the linear propagation of P. 

Implementing the filter with Hill’s model is fairly straight- 
forward. Since Hill’s model is in the relative frame, the shut- 
tle’s position is implicitly included in the states, any param- 
eters derived from the position, can be assumed to be 15 
provided by the GPS receiver’s navigation solution. 

The equation describing the reflected Pseudorange Measure- 
ment is given by: 

20 

h(X)=p r =r sh +r hg -r gs 

The subscripts stand as the Shuttle to Hubble, Hubble to a 
GPS satellite, and A GPS satellite to the Shuttle respectively. 
The GPS satellites position is known from the ephemeredes, 
and Shuttle’s position is known from the receiver’s naviga- 25 
tion solution. If we define, A r =X s -X h , that is the relative 
position between the Hubble and Shuttle and rotate into Hill’ s 
frame we arrive at: 

p=\A r \+\A-X g \-\X g \ 30 

A r is our relative position vector to Hubble in the ECI frame, 
and part of the state vector. Now taking partials with respect 
to A r we arrive at: 


dp r A r Ar — Xg 
3A; = jA^ + |A r -X s | 
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Taking the partials with respect to the state Ar: 


A2=a v .fi 

dA r A r 


A? 

IAT 


+ (Av-V,) 




X g \ 


+ 


±r~X g ) 

|A r -X,| 3 J 


Now, taking the partials with respect to the state Av: 


dp _ Ar A r — Xg 

dT v = \KA~\A r -x g \ 


Now to form the H vector: 


Hp r {x) = 


dp dp 
~dA r dA~ v _ 


For the nonlinear full dynamics implementation, measure- 
ments are given by the measurement equations as a function 
of the states of the system, such that y=h(x). In this formula- 
tion, there are 4 different measurements, pseudorange, differ- 
enced reflected pseudorange, pseudorange rate, and reflected 
pseudorange rate. In order to use these measurements in our 
Kalman filter design, we must linearize the measurements in 
order to form the canonical H matrix. 

The direct psuedorange measurement is given by the equa- 
tion: 


h p (x)=p=r+c(t u +bt+bt D ) 

neglecting the ionospheric, receiver, ephemeris and other 
errors (6+6T rf ), the measurement can expressed in terms of a 
range and bias. 


h p {x) =p = r + b = (x s - x g ) 2 + (y s - y g ) 2 + (Zs~ Z g ) 2 + b 


This forms our measurement vector: 


H(x) = 


dpr 

dA rlx3 


01x3 


where b is the receivers time offset (in meters) from the GPS 
45 constellation. 


The equation describing the reflected pseudorange rate 50 
measurement is given by: 


dp 

H pW =/ x 


The partial with respect to the Shuttles position: 


Now, substituting in the line of sight vectors: 


P r = (V s -V h )- 


Xs-X h 

\Xs-X h \ 


+ (V k -V g )- 


X h -X g 

\x h -x g \ 


55 = x s- x g 

dx s (x s - x g ) 2 + (y s - y g ) 2 + (Zs~ Z g ) 2 

let r sg = ^ (x s -x g ) 2 + (y s - y g ) 2 +(Zs~ Z g ) z 


Now substituting in A=X s -X h , and A v =V s -V h : 


Now taking it with respect to each position vector: 


P = A v - 


A, 

Wa 


+ (A v - V g ) ■ 


\^r-Xg\ 


ys -y g zs-z g - 


65 
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The partial with respect to the Shuttle’s velocity: 


Now extrapolating for each of the shuttle position compo- 
nents, we get a 1x3 vector 


The partial with respect to receiver bias: 


- X s Xfj X s Xg X s Xfj X s Xg X s Xfo X s Xg 
f sh Fgs ^ sh f gs f sh F gs 


The partial with respect to the Shuttle’s velocity: 


The partial with respect to Hubble’s position: 


p=o lx3 

dx v 


The partial with respect to the Shuttle’s bias: 


The partial with respect to Hubble’s velocity: 


20 U J!L =0 
dx b 


The partial with respect to Hubble’s position: 


Now forming a H vector for the direct measurement 


H p (x) = [p- 0 U3 0i ri 1 


(x s -x h f + 


(x h -Xg) 2 + 


(y s - yh ) 2 + (Zs- Zh ) 2 y iyh - y g ) 2 + (z h - z g f 


_ , _ , , . , . _ Now substituting in the range definitions, and taking par- 

ihe dillerenced reflected pseudorange provides lnlorma- . , rrmrr , , . , 

+• u i» i + u 4 +i i t i uu •+• ++ tials oi HST s y and z components we arrive at the H vector: 

tion about the absolute shuttle and Hubble position states. Its J t 

measurement equation is given as follows: 35 

k(x)=p =r sh+ r hy -r sy _ x h -x g _ 

The bias terms difference out, as well as many of the other r h g 

common mode errors greatly simplifying the calculation and dp r _ y*-yh yh~y g _ 

the algebra. The subscripts sh, hg, and gs signify the range 40 dx h r sh r hg 

from shuttle to Hubble, Hubble to GPS transmitter, and Zs-Zh Zh-z g 

shuttle to GPS transmitter. These are equivalently: l r ^g 

r sh=^ {X^ h ) 2 +(y-y h ) 2 +(z-z h f 

45 The partial with respect to Hubble’s velocity: 

r h g =^ (^h-Xg) 2 +Th-y g ) 2 +(^h-Zg ) 2 

V*i fe-Xg^+O.-y^+fz.-Zg ) 2 dr h 

d — = ^ lx3 

In order to form the H vector equation we must first linear- 50 pr 

ize the measurement equation. The partial with respect to the 

Shuttle’s position: XT ~ . TT ^ 

Now forming a H vector for the direct measurement 


Xs-Xf j X s -X g 

(x s - X h ) 2 + I ( X s -Xg ) 2 + 

(y s - yh f + (Zs - Zh ) 2 v (y s - y g ) 2 + (z s - z g f 


Hpr(x) ~ [ „ y P 0^3 — - — 0^3 0 

L °*s\x 3 0*hlx3 


substituting in the definitions of the range vectors we get: 


The equation that relates the Pseudorange rate measure- 
ment to the states is given by: 


dp _ x s -x h _ x s -x g 
dx s r sh r gs 


h(x)=p = (V u -V g )- 


65 
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The measurement is simply the relative velocity along the Because there is no information contained in the pseudor- 

line of sight vector between the GPS transmitter and the ange rate about the state vector of the Hubble, the partials with 

Shuttle. respect to Hubble’s states are 0. 


dp_ 
' dX 


This partials with respect to the scalar x position is given as: 


. to 


dp 

~dX h = 0l * 3 
dp 
dV h 


= 0 ., 


dp . 

£ s =^- v m 


V ( x s - x h ) 2 + (y s ~ yh ) 2 + (zs- z h ) 2 

(x s -x g ) 2 


(V (x s - x h ) 2 + (y s - y h ) 2 + (Zs ~ Zh) 2 ) 

Now substituting in the range formulas to clean up 


dp 

= (V s - V s ) x 


1 (x s -x g y 


r sh r\ h 


Now constructing the H vector for the psuedorange rate 
measurement, we get: 


15 


dp dp 

dXsix 3 dV s 1*3 


20 The reflected pseudorange measurement is essentially the 
sum of the Doppler shift along the paths of the signal. This is 
given by: 

h(X)= p=(V s - V h yi6s sh+ (V h - v g )-[6s hg 
25 

Now, expressing the line of sight vectors in terms of their 
components: 


The same holds for the y and z vectors as well. 


x s - x h x h - x h 

30 J 1 , + (V h -V g )■ 


- X h \ 


\x h -x g \ 


dp 

W S = (V ‘- V ^ 


dp 


1 (y s -r g f 


r sh rl h 

1 ( Zs-Zg ) 2 


rlh 


From this, we need to take the partials with respect to the 
state vector components. For the shuttle position states, in the 
35 x direction: 


Now taking the partials with respect to range, it simply 
leaves the line of sight vectors. 40 


dp [ 1 {x s -x h ) 2 

dx s ys h)x [\x s -x h \ |x s -x A | 3 


Substituting in the range vectors to simplify the equation 


dp 

X s -x h 


we get 





dV s ' 

\X S -X h \ 

45 









dp _ 

= (Vs- 

Vk)x\ 

[J_ 

(x s -x h ) 2 } 

breaking these into components and substituting in the range 


dx s 


4, \ 

equation, we get: 









50 

For the shuttle position states, in 

dp 

x s -x h 







dV^ 

r s h 


dp 

= (Vs- 

V h ) y 

rj_ 

Cys-y/z) 2 ] 

dp 

_y s -yh 

55 

dy s ' 


4s 1 

dVsy 

r s h 









dp _ 

= (Vs- 

V A ) Z [ 

' 1 

te - Zh) 2 1 

dp 

Zs-Zh 


dz s " 

r s h 

ds \ 

dV sl 

r s h 








Assuming a constant receiver bias: 


60 Now taking the partial with respect to the Shuttle’ s veloc- 
ity: 


dp 

dx b 


= 0 


dp X s -X h 


65 
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This simply leaves the line of sight vector to the Hubble; 
expanding this into the x,y, and z components: 


dp = y s -yh 
dV sv r sh 


= -(V s -Vh) x - 


1 (x s - x h f " 


1 (Xfj Xg ) 

\Xh-X g \ ~ |X„-;y 3 


The measurement matrix corrects the Kalman Filter’s 
dynamic estimates with the current measurement estimates. 
The measurement matrix is stacked up in the same order as 
the measurements. 

5 The measurement for the relative filter is given by: 

H=[H pr {X)fl pr {X)] T 

The matrix for the Absolute dynamic filter is given as: 

1Q H=[H p {X)fl p {X)fl pr {X)fl pr {X)] T 

The process noise matrix, Q describes the covariance error 
associated with the model. It generally describes the model 
errors and un-modeled force inputs. In our case, since Q is 
difficult to obtain, it is treated as a tuning matrix, to adjust the 

15 performance of the filter. 


substituting in the range vectors to clean up the equation, and 
extrapolating to the y,z vectors we arrive at: 


1 

(x s - x h ) 2 1 

' 1 

= -(Vs-V h ) x -\ — 

- 3 \ + (Vh-V g ) x - 

— 

[r s h 

d J 

r hg 

... ... ri 


' 1 

M 

y 

i 

g 

T 

- 3 \ + (Vh- Vg) v 



" [fsh 

d J y 

r hg 

... ... n 

(2$ Zfi )1 

' 1 

M 

i 

g 

T 

5 — + v*-v,V 

— 


d J z 

r hg 


Now taking the partial with respect to the Hubble’s veloc- 
ity, we simply get a relation involving the line of sight vectors. 


dp X s -X h X h -X g 
W h = ~ |X S -X*| + 


expanding it into the x,y,z components and substituting in the 
range definitions we get: 


o p is the standard deviation in Hill’s relative position model, 
o v is the standard deviation in Hill’s relative velocity model. 


0 cr^hhxs 

0 0 

0 0 


0 othhxs 0 

0 0 d 


o ps is the standard deviation in the Shuttle’s position model, 
o vs is the standard deviation in the Shuttle’s velocity model. 
o ph and o vh represent the standard deviation in Hubble’ s posi- 
tion and velocity models respectively. 

The R matrix represents the variance in the measurements. 
Ideally the variance is uncorrelated white gaussian noise. 
Since we can estimate the variance of the measurements with 
relative accuracy, this matrix is essentially fixed. 


dp _ x s -x h + x h -x g 
d Vhx r sh r hg 

dp _ y s -yh + yh-y g 

d Vhy r sh r hg 

dp _ Z s ~ Zh + Zh ~ Zg 

d Vhz r sh r hg 


50 o D and o D9 are the differenced psuedorange and reflected 
pseudorange rate measurements. 


Since the reflected pseudorange rate does not involve the 
bias, it is zero. 


Forming the linearized measurement vector we get: 


r dpr dpr dpr dpr 
^ ^ _dX slx3 d V six3 d X hlx3 dV hlx3 _ 


o p is the standard deviation in the pseudorange measurement, 
cfp is the standard deviation in the pseudorange rate measure- 
ment. o D p and o D p are the differenced psuedorange and 
65 reflected pseudorange rate measurements. 

Implementation of a preferred embodiment will now be 
described. 
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The use of reflected GPS signals as a navigational tool may 
be particularly useful when navigating toward an orbiting 
target such as the Hubble Space Telescope (HST) in order to 
perform repairs, replace critical parts or otherwise service an 
orbiting device. The ability to passively and autonomously 5 
navigate to perform servicing missions could lead to 
unmanned service missions as well as reduced weight and 
power consumption. The following description relates to 
navigating the space shuttle to the HST. 

The space shuttle will incorporate a GPS receiver such as to 
that disclosed and describe in US 2006-0082496 to Winter- 
nitz et al. which is incorporated herein in its entirely. The ’496 
receiver is particularly suited for implementing the method of 
the present invention using reflected GPS signals by having 
the ability to acquire and track weak GPS signals. The navi- 15 
gation system will include two antennas; a Right Handed and 
Left Handed Circular Polarized (LHCP, RHCP) antennas. 
The GPS receiver will be powered on as it approaches the 
HST within 1000 m from Hubble. The Shuttle’s position can 
be estimated with conventional methods of navigation using 20 
direct GPS signals and a dynamic orbit model. The receiver 
onboard the Space Shuttle will be acquiring the direct signals 
as it approaches the HST which will be received by the RHCP 
GPS antenna. A convenient feature about a RHCP signal is 
that, when reflected off of a surface such as the HST, it 25 
reverses polarization and becomes left hand circularly polar- 
ized. This gives it a nearly perfect rejection for a reflected 
signal. Therefore, the LHCP Antenna is well suited to pick up 
the reflected signals. 

Both of the antennae will receive the raw RF stream mixed 30 
down to IF, 35 .42 MHz, and then sampled at 8. 1 92 MS/s. This 
data will be processed as previously described to compare the 
direct and reflected signals to estimate the position of the 
Shuttle relative to the HST. 

The use of ever increasing speeds of computers, micropro- 35 
cessor and FPGAs, complex signal processing algorithms 
will be used to apply the aforementioned filtering techniques 
and signal acquisition methods to determine the structure and 
composition of the reflected signals. 

The GPS receiver used to record and process the data is the 40 
Navigator GPS. The GPS receiver is running at 65.536 MHz 
using a General Dynamics Coldfire processor and has three 
application specific Actel AX2000 FPGA’s. The enhanced 
capabilities of this receiver are aided by the application spe- 
cific FPGAs that implement an algorithm developed to 45 
acquire weak signals in real time as disclosed ’496 to Win- 
ternitz et al. This algorithm is an FFT based acquisition algo- 
rithm that operates in the frequency domain. It is able to 
acquire and track signals down to 25 dB-Hz. This is approxi- 
mately 1 0 dB more sensitive than current off the shelf GPS 50 
receivers. This will enable the Navigator to acquire the 
weaker reflected signals off of the HST. 

The Navigator receiver will be configured to track 24 chan- 
nels, 12 of them connected to the LHCP antenna and the 
remaining 1 2 connected to the RHCP antenna. While tracking 55 
a direct signal, the reflected channels will be directed to 
search in a time delayed space with approximately the same 
Doppler frequency. This will tremendously aid acquisition of 
the weaker reflected signals. 

The aforementioned extended Kalman filter will run on the 60 
navigation processor taking in the RHCP and LHCP measure- 
ments. Thus the Navigation system will be able to estimate 
the position of the Shuttle in orbit and the relative position of 
the HST by comparing the directed and reflected signals of a 
plurality of GPS satellites as previously discussed. Using a 65 
receiver to receive both direct and reflected GPS signals and 
comparing/processing the embedded ranging data (time 
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delay etc.), coupled with the dynamic model in the extended 
Kalman filter, the navigation system of the present invention 
can passively and autonomously track a relative position of an 
orbiting spacecraft (Shuttle) to an orbiting target (HST). The 
ability to passively and autonomously passively navigate a 
space craft toward a target not only reduces weight and power 
consumption, but also paves the way for potential precise 
passive navigation of space craft for unmanned service mis- 
sions. 

A simulation and implementation of the method of the 
present invention will now be discussed in regards to a Shuttle 
mission approaching the Hubble Telescope. 

The Hubble Rendezvous truth- scenarios are chosen such 
that the initial conditions will result in a near collision of the 
Shuttle and Hubble. The truth model is generated with STK 
using a high fidelity orbit propagator. 

The Space Shuttle may approach Hubble from the aft end 
with its bay doors pointed towards Hubble. This allows the 
Shuttle to dock cleanly with Hubble. 

During the servicing mission, the left and right handed 
GPS Antennae will be placed on the MULE inside the bay of 
the shuttle looking outwards. 

The model focuses on simulating the uncontrolled motion 
of the two bodies. In order to represent the “docking”, initial 
conditions were chosen such that the two space craft come 
within close proximity to each other. 

The model was created using STK, given a set of initial 
conditions, and propagated over an orbit period. The STK 
model chosen was a 2 1 -order gravity model with atmospheric 
drag on both the Shuttle and Hubble, the coefficients of drag 
chosen as to represent the relative mass and area ratios. 


Units 

Shuttle 

ECI Position 

(-4673.7, -4525.4, 2443.5) T 

Km 


ECI Velocity 

(5.604, -4.472, 2.43 8) T 

Km/s 


Mass 

109,000 

kg 


Area 

58.4 

m 2 


Q* 

1.2 


Hubble 

ECI Position 

(-4673.9, -4525.9, 2443.7f 

Km 


ECI Velocity 

(-5.604, -4.471, 2.437)'* 

Km/s 


Mass 

11,110 

kg 


Area 

15.4 

m 2 



0.8 



The area for the Shuttle was approximated using the circular 
cross sectional area of the Shuttle. The same was done for the 
Hubble. 

The visibility of each GPS satellite is determined by the 
geometry of the GPS constellation with respect to the location 
of Hubble and the Shuttle. The satellite visibility is essential 
to determining how many measurements both direct and 
reflected are available. This number of measurements avail- 
able influences the absolute and relative position accuracies. 

The GPS constellation is constructed and propagated using 
the keplerian elements given by an almanac file. The almanac 
file is generated and uploaded to each GPS satellite. It repre- 
sents a limited set of the keplerian orbits, in order to provide 
a rough position estimate of the all of GPS satellites. These 
orbits are then propagated in time, providing a good estimate 
of the satellites orbits. This is purely used to generate a mov- 
ing rough estimate of the GPS constellation for the simula- 
tion. 

In practice the GPS constellation is closely monitored, and 
each satellite has a new ephemeris updated every 3 hours. In 
addition to updating their GPS Almanac files every 1 2 hours. 
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The Almanac parameters are given by: 
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plate. In physical optics, the equation describing the RCS for 
a metallic circular disk is: 


Parameter 

Symbol 

Scaling 

Time of Almanac 

toa 

secs 

Eccentric Anamoly 

E 

none 

Inclination 

iO 

semicircles + 0.94247 

Rate of Right Ascention 

Q 

semicircles 

Semi-major axis 

AV 2 

none 

longitude of perigee 

X 

semicircles 

Argument of Perigee 

(0 

semicircles 

Mean Anamoly 

M 

semicircles 


Since we are only interested in observing those GPS Sat- 
ellites with a reflected pair, and since the relative distance is so 
small, we can assume that any GPS satellite seen by the space 
Shuttle can also be seen by the Hubble. 

The visible signals are done by masking out those GPS 
satellites whose main lobes cross through the earth, or are not 20 
in view of the shuttle. This is done by calculating the angle 
between the line of sight vector, and position vector of the 
shuttle as seen in FIG. 4 . If (3 is greater than 13.9° and less than 
2 1 .3° then it is within the main lobe andjust over the lim of the 
earth. If |3 is less than 13.9°, then we need to check to make 25 
sure the signal does not pass through the earth. Thus, the dot 
product of the two vectors should be negative. 

The bi-static radar equation described in is given as: 


R r 

R* 

P r 

G r 

G* 

F r 
F* 

K 

T„ 

D„ 
(S/N) 
L r 
L * 


P7 -GtG b aV b f|-f| y a 
{An ) 2 KT.s B jy ( S/N) mj - n L7L /? y 


Transmitter to target range 

Target to receiver range 

Transmitter power 

Transmit antenna power gain 

Receive antenna power gain 

Bistatic radar cross section 

Pattern propagation factor for transmit to target 

Pattern propagation factor for receiver to target 

Boltzman’s constant 

Receiver noise temperature 

Noise bandwidth of the receivers pre-detection filter 
signal-to-noise ratio required for detection 
Transmit system losses 
Receive system losses 
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40 


45 


This equation is used to determine the maximum range to 
the target. However, rearranging terms, it is possible to pose 
this in the form of SNR received. 

There are several universal constants from this particular 
implementation. X is the wavelength of the GPS LI signal. 
The GPS LI signal is transmitted at 1575.42 MHz, giving it a 
wavelength of 0.1903 m. Boltzmann’s constant is given by 55 
$1.3806503xl0 -23 J/K. F r and F^ are the pattern propagation 
factors, these factors describe the loss due to transmission in 
an unclean environment. F r and F^ can be estimated at around 
2 dB total for the ionospheric effects. 

The bi-static radar cross section is calculated in much the 60 
same way as the standard radar cross section. However, the 
transmit angle needs to be taken into account, as well as 
properties limited only to the bi-static scenario. 

Since the Hubble can be approximated as a cylinder and the 
approach profile is towards the aft-end. The Hubble can 65 
approximated by the reflective properties of a circular flat 
plate combined with a cylindrical surface normal to the flat 


Acos0 

1 f\7T 

Ji{kdsmQ) 2 

A 

kdsinO 


Where A is given as the physical area of the disk, d is the 
diameter, and J : (x) is the Bessel function of the first kind of 
order one. X is the wavelength. 

The Hubble space telescope has a diameter of roughly 4.2 
m. Thus giving an area of roughly 13.85 m 2 . 

For normal incidence, this equation simplifies to: 


CT = 


An A 2 
~AA 


This would represent the monostatic cross sectional area of 
the Hubble space telescope with the observer looking perpen- 
dicular to the aft end. Thus, the maximum possible cross- 
sectional area, looking at the aft end is 36.81 dB-sqm. It is 
interesting to note that the GPS signal will get gain from the 
reflections. If the observer was oriented perpendicular to 
Hubble, it would appear as a cylinder to the viewer. In this 
case, the RCS is given by: 


2 sm(klsm(theta )) 2 

cr = kal — . — 

klsm{theta) 


Where a=2.1 m is the radius, 1=14.2 m is the length of the 
Hubble, k is the wave number, and X is the wavelength. This 
provides a predictable radar cross section. With the assump- 
tion since at the extremes of both of these figures, the RCS 
drops off significantly, and effectively take the maximum 
RCS. Of the two plots, the side of the cylinder offset by 90° to 
form a combined RCS model for the cylinder over 1 80° . Since 
a cylinder is symmetrical, this is effectively models the full 
cylinder over 360°. 

In order to account for the bi-static nature of the problem, 
it is necessary to define a new angle called the bi-static angle. 
This angle is defined as the difference between the receive and 
transmit antennae, p=0 r -0„ the bisection of this angle is used 
to approximate the Bi-static RCS, o b . The bi-static RCS can 
be grouped into 3 sections: the pseudo -monostatic RCS 
region, the forward scatter RCS region, and the bi -static RCS 
region. In the pseudo -monostatic region, the Crispin and Sei- 
gal mono static -bi -static equivalence theorem predicts that for 
a sufficiently small angle, the RCS is equivalent to the mono- 
static RCS measured the bisector of the bi static angle. This 
theorem applies to a sufficiently smooth surface, such as a 
sphere, a cylinder, or cone. Some empirical data for a cylin- 
der, performed at 35 GHz, the pseudo-monostatic region 
extends to roughly 20°, while the bistatic is within 
20°^(3^140° and forward scatter at (3>140 °. 

The bi-static region begins where the pseudo-monostatic 
region ends, it is where the theorem fails to predict an accurate 
RCS. Unfortunately, there is no empirical formula for this 
equation. Experimental data shows a downward trend, from 
-2 dB to -15 dB, in this case. It should suffice to use a linear 
approximation (in dB). 

The third bi-static region is the forward scattered region, 
this occurs when the bi-static angle approaches $ 1 80"o$. The 
forward scatter can be approximated by treating the shadow 
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area as a uniformly illuminated antenna aperture. The roll off 
can be approximated when fl-|3 is substituted for the angle off 
aperture normal. This function closely matches J 0 (x)/x where 
J 0 (x) is the Bessel function of order 0. In our situation, many 
of the GPS satellites will lie in this area, due to being below 5 
the constellation. This is to our benefit, allowing us to pick up 
signals, being at such an odd orientation. Luckily due to the 
change in polarization of the reflected signal, this will provide 
us with isolation from the direct signal from the satellite. 

From these bi -static regions we can now form a full to 
approximation for the bi-static radar coefficient. The plot 
below represents a b with the flat face of the cylinder oriented 
towards the receiver. This should provide a reasonable fit for 
a perfectly reflecting disk. 

The antennae chosen for this mission are two hemispheri- 15 
cal antennas. The antenna gain is given as 3 dBic. Built into 
each antenna is the low noise amplifier. The LNA has a gain 
of 26 dB with a noise figure of 2.8. This can be used to 
calculate the system noise temperature. This equation is given 
as: 20 

T=10 log(T A (NF-l)) 

Where, T 4 is the temperature at the antenna. Since the LNA 
is at the antenna, the line losses are included in the analysis. 
From the spec sheet, the NF of the antenna is given as 2.8 dB. 25 
A conservative approximation for this noise figure is to use 
the temperature at the antenna, T A , as room temperature or 
290 K. This approximation holds because the antennae are 
lying on the MULE. The MULE will be heated in order to 
keep the electronics from freezing. From this equation we 30 
find that the receiver noise temperature is 24.4 dB-K. 

The L^ and L r and receiver and transmit system losses 
respectively will be taken into account as a total system loss. 

A general approximation of loss in a communications system 
is about 2 dB. Usually this value is measured quantitatively, 35 
but, without a complete system to test, it is difficult to arrive 
at a true figure. The noise bandwidth, B n is generally taken as 
1/T, where T is the pre-integration period. For Navigator, it 
has an adjustable preintegration period, however the nominal 
case is 1 ms. Other losses occur such as Polarization mis- 40 
match contributing 3 dB and in the case of a reflection, there 
is an efficiency which is dependent on the reflectivity of the 
object. In the case of the Hubble, which is coated in mylar, it 
will have a high reflectivity coefficient, this is typically 
90-98%, depending on the quality. 45 

The range calculations R r is simply the geometric distance 
between the Shuttle and the Hubble, while R r is simply the 
geometric distance between each GPS Satellite and the 
Hubble. Factoring the free space loss due to range the equa- 
tions: 50 


R T (dB) = 201o 


/?K(dB) = 201o 
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At each time-step, for each visible GPS satellite, each 
range value is calculated. With the assumption that the Shuttle 60 
is oriented with its bay doors towards the aft bulkhead of 
Hubble, it is possible to estimate the angle p and which 
surface of the cylinder the reflection is striking. The truth 
positions were generated from a file that is the simple 2-body 
relative motion, coupled with some thruster firing for a real- 65 
istic closing profile. This produces an expected satellite vis- 
ibility given the aforementioned conditions. It should be 


noted that the model does not take into account the reflections 
from the solar panels and assumes an orientation with the 
antennae directly facing the aft bulkhead. 

The direct signal is generated using the truth positions of 
the Shuttle and Hubble. A receiver bias is added on top of the 
signal, as well as the ionospheric error given by the equation 


S2.1TEC 

®iono — i = 

sin 2 (EL) + 0.76 + sin 2 (EL) 

TEC is the total electron content, given by 2el7. F c is the 
carrier frequency 1 .57542 GHz, and EL is the elevation of the 
GPS satellite over the horizon. 

The direct signal pseudorange rate is generated by the rate 
of change of the line of sight vectors. In practice, most GPS 
receivers rely only on the pseudorange measurements or use 
carrier phase to generate the velocities. The pseudorange rate 
measurements are generally noisy and actually degrade the 
performance of the filter. 

A variance of 5 meters was chosen for the direct pseudo- 
ranges, and 25 cm/s for the pseudo -range rates. It is composed 
of the clock drift, ephemeris errors, phase noise, and other 
unmodeled sources. The ionospheric bias was neglected in 
this simulation, as it would have been differenced out in the 
measurement generation. The pseudo-range variance may be 
increased by a nominal amount, 0.25 m, to reflect the iono- 
spheric model discrepancies. 

The differenced reflected measurements were constructed 
from the total path lengths between the GPS Satellite, 
Hubble, and the Shuttle. The noise added is from the loss in 
signal quality due to the reflective properties and small iono- 
spheric error. Since the measurement is a differenced 
reflected signal with that of the direct path, the bias and the 
ionospheric terms become small and are negligible. In addi- 
tion to the ionospheric error practically canceling out, the 
front end bias, clock bias, thermal jitter and other common 
mode errors also disappear. 

A variance of 25 meters was chosen for the reflected 
pseudo -ranges, and 50 cm/s for pseudo -range rates. Although 
the true variance is unknown, it is extremely difficult to model 
on the earth. Future tests will include an anechoic chamber 
testing or data reduction of the rendezvous experiment. 

While the method and apparatus of the present invention 
has yet to be tested oil a service mission, simulations have 
proven the viability and utility of the use of comparing 
reflected and directly received GPS signals to passively navi- 
gate spacecraft to an orbiting target. The present invention is 
currently scheduled to fly in August 2008 on Hubble Servic- 
ing Mission 4. A further discussion of such simulations can be 
found in Ian Cohen’s thesis entitled “Relative Navigation for 
Hubble Mission Using Reflected GPS signals”, published by 
the University of Maryland in May 2007 and is hereby incor- 
porated herein by reference in its entirety. 

One additional use of the method of using reflected GPS 
signals as a passive navigation system is to compliment exist- 
ing navigation tools. This passive system can not only provide 
an extra measurement tool offering additional redundancy, it 
could also be used to navigate a space craft to within sufficient 
proximity of an orbiting target reducing the amount of time 
active systems such as LIDAR and RADAR are powered. 
Because this method and system of the present invention may 
be employed by simply adding the second LHCP antenna and 
loading the algorithmic complexity onto existing systems, 
redundancy and reduced power consumption can be achieved 
with little modifications to existing equipment. 
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The foregoing description of the specific embodiments will 
so fully reveal the general nature of the invention that others 
can, by applying current knowledge, readily modify and/or 
adapt for various applications such specific embodiments 
without departing from the generic concept, and, therefore, 5 
such adaptations and modifications should and are intended 
to be comprehended within the meaning and range of equiva- 
lents of the disclosed embodiments. It is to be understood that 
the phraseology or terminology employed herein is for the 
purpose of description and not of limitation. Therefore, while 10 
the embodiments of the invention have been described in 
terms of preferred embodiments, those skilled in the art will 
recognize that the embodiments of the invention can be prac- 
ticed with modification within the spirit and scope of the 
appended claims. 15 

The invention claimed is: 

1. A method of passively navigating an orbiting spacecraft 

towards an orbiting target, said method comprising the steps 
of: 20 

receiving a plurality of direct signals, each of said plurality 
of direct signals transmitted from a corresponding one of 
a plurality of celestial GPS earth orbiting satellites trans- 
mitting ranging data signals; 

receiving a reflected signal from said corresponding one of 25 
said plurality of GPS satellites reflected off said orbiting 
target; 

comparing said direct signals and said reflected signals; 
and 

determining a position of said orbiting spacecraft relative 30 
to said orbiting target based on said step of comparing 
said direct signals and said reflected signals. 

2. The method of claims 1, wherein said step of comparing 
said direct and reflected signals includes the step of determin- 
ing a time delay between receiving said direct signal and 35 
receiving said reflected signal from a particular GPS satellite. 

3. The method according to claim 1 wherein said steps of 
receiving said direct and reflected and said step of comparing 
said direct signals and said reflected signals are autono- 
mously and passively performed on said orbiting space craft. 40 

4. The method of claim 1, step of determining a relative 

position of said orbiting spacecraft includes tracking the posi- 
tion of the spacecraft and updating the position of the space- 
craft as a function of the tracked position and said ranging 
data. 45 

5. The method of claim 4, wherein step of determining said 

relative position of said orbiting spacecraft includes filtering 
said ranging data in said reflected signals and said direct 
signals. 5Q 

6 . The method of claim 5, wherein step of determining said 
relative position of said space craft includes employing an 
extended Kalman filter to compute the position of the orbiting 
spacecraft relative to said orbiting target. 

7. The method according to claim 1, wherein said step of 55 
receiving said direct signals is performed with a left handed 
circular polarized antenna and said step of receiving said 
reflected signal is performed with a right handed circular 
polarized antenna. 

8 . The method according to claim 1, wherein said step of 60 
comparing said direct and reflected signals includes process- 
ing said signals into pseudorange measurements. 

9. The method of claim 8 , wherein said step of processing 
includes processing said direct signal into a pseudorange 
measurement representing a distance directly between said 65 
orbiting space craft and a corresponding satellite and process- 
ing said reflected signals into a pseudorange measurement 
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representative of a sum distance between said orbiting space 
craft to said orbiting target and said orbiting target to said 
corresponding satellite. 

10. The method of claim 3, further comprising a step of 
implementing a extended Kalman filter to process said rang- 
ing data to generate and update estimated position data 
thereby autonomously generate real time estimates of said 
relative position of said orbiting space craft relative to said 
orbiting target. 

11. A navigation system to passively navigate an orbiting 
spacecraft to an orbiting target, said navigation system com- 
prising: 

a receiver platform adapted to receive a plurality of signals 
carrying ranging data transmitted by at least one of a 
constellation of satellites; said receiving platform 
including; 

a first antenna adapted to receive at least one of said 
plurality of signals direct from a corresponding satel- 
lite; 

a second antenna adapted to receive said at least one of 
said plurality of signals from said corresponding sat- 
ellite reflected off said orbiting target; 

a processor adapted to process said ranging data received 
by said first and second antennas to thereby autono- 
mously and passively locate a relative position of said 
orbiting spacecraft relative to said orbiting taiget. 

12. The navigation system according to claim 1 further 
comprising; 

an extended Kalman filter adapted to process said ranging 
data received by said receiver platform and generate and 
update estimated position data to autonomously gener- 
ate estimates of said relative position of said orbiting 
space craft relative to said orbiting target. 

13. The navigation system according to claim 1, wherein 
said processor is adapted to process said at least one of said 
plurality of signals direct from a corresponding satellite into 
a pseudorange measurement representing a distance directly 
between said orbiting space craft and said corresponding 
satellite and processing said at least one of said plurality of 
signals reflected of said orbiting target into a pseudorange 
measurement representative of a sum distance between said 
orbiting space craft to said orbiting taiget and said orbiting 
target to said corresponding satellite. 

14. The passive navigation system according to claim 1, 
wherein said first antenna is a left handed circular polarized 
antenna and said second antenna is a right handed circular 
polarized antenna. 

15. The passive navigation system according to claim 4, 
wherein each of said first and second antennas is a hemi- 
spherical antenna each incorporating a low noise amplifier. 

16. The passive Navigation system according to claim 13, 
wherein said receiving platform receives signals transmitted 
from at least four satellites and said processor implements 
said extended Kalman filter to measure a relative position of 
said orbiting spacecraft relative to said orbiting target. 

17. A method of passively navigating a moving body 
towards a moving target, said method comprising the steps of: 

receiving a plurality of direct signals, each of said plurality 
of direct signals transmitted from a corresponding one of 
a plurality of celestial GPS earth orbiting satellites trans- 
mitting ranging data signals; 

receiving a reflected signal from said corresponding one of 
said plurality of GPS satellites reflected off said moving 
target; 

implementing an extended Kalman filter to process said 
ranging data to generate and update estimated position 
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data thereby autonomously generate real time estimates 
of said relative position of said moving space craft rela- 
tive to said moving target; 

comparing said direct signals and said reflected signals; 
and 
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determining a position of said moving body relative to said 
moving target based on said step of comparing said 
direct signals and said reflected signals. 



